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ABSTRACT 

We propose that inward, subsonic flows arise from the local dissipation of turbulent motions in molecular 
clouds. Such "turbulent cooling flows" may account for recent observations of spatially extended inward motions 
toward dense cores. These pressure-driven flows may arise from various types of turbulence and dissipation 
mechanisms. For the example of MHD waves and turbulence damped by ion-neutral friction, sustained cooling 
flow requires that the outer gas be sufficiently turbulent, that the inner gas have marginal field-neutral coupling, 
and that this coupling decrease sufficiently rapidly with increasing density. These conditions are most likely met 
at the transition between outer regions ionized primarily by UV photons and inner regions ionized primarily by 
cosmic rays. If so, turbulent cooling flows can help form dense cores, with speeds faster than expected for 
ambipolar diffusion. Such motions could reduce the time needed for dense core formation and could precede 
and enhance the motions of star-forming gravitational infall. 

Subject headings: ISM: clouds — MHD — stars: formation — turbulence 



1. INTRODUCTION 
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! Q Recent maps of molecular cloud cores indicate inward mo- 
tions that appear too extended to match models of gravitational 
infall and too fast to match models of gravitational contraction 
limited by ambipolar diffusion. In L1544 (Tafalla et al. 1998) 
and other "starless" cores (Lee, Tafalla, & Myers 1998) and in 
some "turbulent" cores with embedded stars (Mardones et al. 
1998), a common feature is a zone of radius 0.05-0. 1 pc, where 
^^most of the observed spectra have "infall asymmetry" — the 
skewing of an optically thick line to the blue of an optically 
OO thin une m a region with inward motions and centrally elevated 
Q^ excitation temperature (Leung & Brown 1977). This zone ex- 
""^ tends beyond the dense core as traced by the half-maximum 
contour of integrated intensity of the 1-0 line of N 2 H + , sug- 
i gesting that the motions are associated with the formation of 
Othe core. 

+3 For LI 544, a detailed analysis of the spectral profiles (Tafalla 
^5et al. 1998) indicates that the observations do not match the 



properties of "inside-out" gravitational collapse (Shu 1977). In 
that model, a star of luminosity several L should have formed, 
whereas no star more luminous than 0.1 L Q is known. The 
observed speeds, ~0.1 km s~\ do not match the gravitational 
collapse of a uniform sphere from rest (Spitzer 1978), which 
predicts much less condensation than is observed. The inward 
motions may be inconsistent with ambipolar diffusion (Mestel 
& Spitzer 1956; Mouschovias 1976) since they are seen in both 
ionic and neutral lines and since their speed at 0. 1 pc exceeds 
that predicted by the model "UV B" of Ciolek & Mouschovias 
(1995) by a factor of ~5. 

These comparisons do not rule out all versions of gravita- 
tional infall or all versions of gravitational contraction limited 
by ambipolar diffusion (e.g., Safier, McKee, & Stahler 1997; 
Basu 1998; Li 1998). Nonetheless, they suggest that nongrav- 
itational motions should also be considered as ways to form 
dense cores. In this Letter, we propose that some extended 
inward motions reflect pressure-driven cooling flows associated 
with local dissipation of turbulence. In a related paper, we 
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consider additional mechanisms of inward motions (Lazarian 
& Myers 1998). 

2. TURBULENT COOLING FLOWS 

"Turbulent cooling flows" may follow localized dissipation 
of turbulence in molecular clouds, as "cooling flows" follow 
radiative cooling in clusters of galaxies (e.g., Cowie & Binney 
1977). Molecular cloud turbulence likely has a substantial 
MHD component (Arons & Max 1975), since the Alfven speed 
is comparable to the nonthermal component of the velocity 
dispersion (Myers & Goodman 1988; Crutcher 1998). MHD 
waves and turbulence can dissipate due to wave damping by 
ion-neutral friction (Kulsrud & Pearce 1969), nonlinear wave 
steepening (Cohen & Kulsrud 1974; Zweibel & Josafatsson 
1983), and interactions of nonlinear Alfven waves (Kraichnan 
1965). MHD waves excited during core formation are expected 
to dissipate due to ion-neutral friction within a few core free- 
fall times unless new waves are internally excited, e.g., by 
winds from young stars (Nakano 1998). Compressible MHD 
turbulence decays nearly as fast as nonmagnetic hydrodynamic 
turbulence, in little more than a free-fall time, according to 
simulations (Ostriker, Gammie, & Stone 1998; MacLow et al. 
1998). A turbulent cooling instability may arise as density in- 
creases, reducing the field-neutral coupling and velocity dis- 
persion, causing an inward pressure gradient and flow, increas- 
ing the density further (Myers & Khersonsky 1995). Thermal 
cooling in molecular clouds (Meerson, Megged, & Tajima 
1996) is closely related, as are cooling to form clouds 
(Elmegreen 1989) and wave-driven flows in diffuse clouds (El- 
megreen 1997). 

We present a simple model to specify the conditions that 
allow a turbulent cooling flow. We assume a uniform outer 
region of density p out and a uniform spherical inner region of 
radius r and density p in . The regions have common isothermal 
sound speed <j T and common Alfven speed v A , in accord with 
the tendency for magnetic field strength to increase as p" 2 
(Crutcher 1998). The outer and inner pressures on the boundary 
surface normal to the field are respectively P om = p oul (a^ + 
^NT.out) and -fL = Pi„(°r + ^NT.inX where a NT denotes the non- 
thermal component of the velocity dispersion. Each region can 
propagate MHD waves, and the outer region has sufficient stars 
and stellar winds to maintain waves with equipartition energy 
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Fig. 1. — MHD wave power transmission / into a uniform spherical region 
as a function of W, the ratio of radius to MHD cutoff wavelength. The upper 
and lower curves represent monochromatic transmission at the lowest and 
highest frequencies of wave propagation in the region. The middle curves are 
averages over a wave energy spectrum ~oT 1/2 as for the inertial range of 
uncorrelated MHD turbulence, assuming that w mllI is set either by the end of 
the turbulent inertial range or by the critical frequency for ion-neutral damping. 
These averages if) indicate significant damping of MHD waves in dense cores 
(ionized by cosmic rays) and in their surrounding molecular clouds (ionized 
primarily by UV photons). 



proper average over all radii and angles, and found that / lies 
within 15% of / for < u < 3. 

So far, this model applies to any kind of turbulent motions 
that exert pressure and to any dissipation mechanism. To make 
specific predictions, we consider the example of dissipation 
due to wave damping by ion-neutral friction, since this mode 
of turbulent dissipation appears likely to dominate over shocks 
and nonlinear mode coupling in the relatively quiescent regions 
in which extended inward motions are currently observed. We 
evaluate k from the dispersion relation for waves of angular 
frequency co in a lightly ionized magnetized gas (Kulsrud & 
Pearce 1969): 



v A [l + (cot) 2 ]> 



sin 



1 



- arctan (cot) 



(3) 



where t is the neutral-ion collision time. The angular frequency 
ranges from co mill « 2Trv A /r, which corresponds to the longest 
wavelength wave that can propagate in the inner region 
(X = r), to the cutoff frequency co max . We use two values that 
probably bracket co max — co ra „(turbulent), corresponding to the 
high-frequency end of the inertial range for MHD turbulence, 
beyond which the spectrum decays exponentially (e.g., Bis- 
kamp 1993), and co max (ion-neutral), where the wavelength 
equals the ion-neutral cutoff wavelength X (Kulsrud & Pearce 
1969). In terms of the field-neutral coupling parameter W = 
r/\ (Myers & Khersonsky 1995), co max (turbulent) « co ram W 2 ' 3 
and co m „(ion-neutral) = co min W. We assume a wave energy spec- 
trum E(o)) do) ~ or 3 ' 2 do) as for the inertial range of uncorre- 
lated MHD turbulence (Kraichnan 1965), so 



J "mi 



do) co 3/2 exp (-2/cr) 



do) co 



(4) 



density, so that a^ T?oat = (v A )/3. The inner region has no internal 
sources of wave excitation, so its turbulence arises solely from 
externally incident waves. The inner turbulence is sufficiently 
"wavelike" that it has negligible shocks and nonlinear inter- 
actions among its waves. A similar model of MHD wave pres- 
sure is presented by Martin, Heyvaerts, & Priest (1997). 

The net inward pressure is AP = P oM — P in , or from the above 
equations 

AP = [p out a 2 ][l " V + (* 2 /3)(l " V8)l (1) 

where rj = pJp out ,x = v A /a T , and g = a 2 T-m /a 2 T-out . If AP > 
0, the resulting inward flow has speed v « [AP/p out ] 1,2 , and con- 
sequently the flow has sonic Mach number 

»/ff r «[l -ij + (* 2 /3)(1 -V8)V' 2 - (2) 

From equation (2), the density contrast r\ reaches its maximum, 
1 + x 2 /3, when v = and ij^<1. 

To obtain the nonthermal velocity variance ratio g, we as- 
sume an isotropic flux of externally incident waves, which have 
negligible reflection at the boundary. We take the spatial mean 
of the monochromatic wave power transmission from the sur- 
face to the inside of the sphere as/ = e~", where u = 2/cr and 
k is the imaginary part of the uniform propagation constant. 
We checked that / is a good enough approximation to /' = 
3iTV"[(l + 3m" 1 ) cosh (k) - (1 + u 1 + 3iT 2 )sinh («)], the 



From equation (3) and the adopted definitions of frequency 
limits, equation (4) can be expressed in terms of the single 
variable W. We evaluated equation (4) by numerical integration 
for W < 10. For larger W we used the approximation cot <C 1, 
which conveniently reduces equation (4) with sufficient ac- 
curacy to a combination of elementary functions and the in- 
complete gamma function. 

Figure 1 shows that g decreases with decreasing field-neutral 
coupling parameter W, depends sensitively on the frequency 
range, and falls below unity even for "marginal" coupling, i.e., 
W< 100. For example, when W = 100, g = 0.9 (negligible 
damping) if only the lowest frequency waves, with co = co min , 
are considered; g = 0.35-0.40 (substantial damping) for the 
spectral averages in equation (4), andg « (complete damping) 
if only the waves at either value of co max are considered. For 
the spectral averages, g decreases from ~1 when log W> 3, as 
for regions of density n < 10 3 cirT 3 (Myers & Khersonsky 
1995), to ~0.1 when W = 10-30, as for dense cores with 
n = 10 4 -10 5 cm" 3 (Williams et al. 1998; Bergin et al. 1998), 
to ~0 when W = 1 as for completely cutoff condensations. 
These results change only slightly if the wave energy spectrum 
varies as oT 5/3 , as expected for hydrodynamic turbulence (Kol- 
mogorov 1941) or "strong Alfvenic turbulence" (Goldreich & 
Sridhar 1995). These results support the idea that starless cores 
with low ionization should soon dissipate the turbulence as- 
sociated with their formation (Nakano 1998). 

Figure 2 shows requirements for "instantaneous" and "sus- 
tained" turbulent cooling flow. For our example of ion-neutral 
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Fig. 2. — Critical conditions for turbulent cooling flow. For each value of 
the density contrast rj between inner and outer regions and of the Alfven Mach 
number x, a turbulent cooling flow can exist only when the field-neutral cou- 
pling parameter W in the inner region is less than the critical curve shown. 
Thus, instantaneous flow requires significant turbulence in the outer region 
(x> 1) and marginal coupling in the inner region (W < 100), and sustained 
flow also requires sufficient decrease in coupling with increasing density. The 
data points indicate mean values of W and x for samples of low-mass and 
more massive cores, and each error bar represents the standard error of the 
mean. These data are consistent with cooling flows onto dense cores. 

damping, the condition for flow at any moment, AP > in 
equation (1), implies that the coupling parameter W be less 
than an upper limit W max given by equations (1) and (4) for 
specified values of x and rj. Figure 2 shows log W max as a 
function of x for density contrast values ij = 1.5, 2, and 3. In 
order to have instantaneous cooling flow into a region with 
substantial density contrast, e.g., rj = 2, W < W max implies that 
the outer region must be turbulent (x >~ 2) and the inner region 
must have "marginal" field-neutral coupling (W < 100). Fur- 
ther, if the flow is to be self-sustaining, W must decrease suf- 
ficiently rapidly as the inner density increases, e.g., for x = 
3, from W< 300 when rj = 1.5 to W< 10 when rj = 3. 

A likely site for such a sustained turbulent cooling flow is 
an outer region dominated by U V ionization and an inner region 
dominated by cosmic-ray ionization, with the boundary at vi- 
sual extinction A vfi « 4 mag from a cloud surface illuminated 
by the standard interstellar radiation field (McKee 1989). This 
is so because the inward pressure in equation (1) is greatest at 
radii ~0.1 pc for low-mass cores, according to the ionization 
model of McKee (1989) and the density model of Caselli & 
Myers (1995), as in Myers (1997). This radius is similar to 
that of observed infall asymmetry. At larger radii, the coupling 
W is probably too strong for significant flow. At smaller radii, 



W is weak enough for flow to arise, but W may not decrease 
fast enough with increasing density to sustain the flow. These 
conjectures need more detailed models for confirmation. 

Figure 2 also indicates that observational estimates of W and 
x are consistent with requirements for turbulent cooling flows 
onto the typical dense core. Each of the two data points rep- 
resents a mean value of W, according to determination of elec- 
tron abundance from observations of H 13 CO + , DCO + , and C ls O 
and from chemical models of neutral and ion abundance, for 
samples of 20 low-mass cores (Williams et al. 1998) and seven 
"massive" cores (Bergin et al. 1998). The value of x is based 
on assuming that the cores have kinetic temperatures 10 and 
15 K, respectively, and that the velocity dispersion in the outer 
region is represented by the mean line width of C ls O, which 
traces lower density gas than do the lines of H 13 CO + and DCO + . 

The mean density contrast t\ = 2-3, which is consistent with 
the models and data in Figure 2, is also comparable to that 
expected for gravitationally bound condensations. The mean 
density over a spherically symmetric region of radius r is 
greater than the density at r by a factor of 2.44 for a critically 
stable Bonnor-Ebert sphere (Bonnor 1956; Ebert 1955) and by 
a factor of 3 for a singular isothermal sphere (Shu 1977). Thus, 
cooling flow may help condense gas to the point at which it 
is self-gravitating and may then accompany the ensuing grav- 
itational infall. If so, cooling flows can modify and enhance 
the process of gravitational infall. 

The characteristic speeds of inflow and boundary expansion 
are similar to the sound speed. We assume x = 3 as in a region 
whose line widths are dominated by turbulent motions and that 
y] has reached half of its maximum value, i.e., r] = 2. Then, 
for W = 20 as in many dense cores, equations (2) and (4) give 
v = l.Ooj, assuming that co max = w(turbulent). Similarly, the 
outward speed of the boundary between the inner and outer 
regions drldt is estimated from the continuity equation if 77 is 
constant as assumed above: drldt = v/rj, or for x = 3, 77 = 2, 
and W = 20, drldt = 0.5o> For 10 K gas, the condensation 
radius would grow from 0.05 to 0.10 pc in ~0.5 Myr. 

3. DISCUSSION 

The pressure-driven flow described here may be an effective 
way to condense molecular cloud gas. Such flows may occur 
wherever turbulent motions dissipate locally, with flow speed 
less than the effective sound speed of the outer gas. For wav- 
elike MHD turbulence with ion-neutral damping and a a/ 3 ' 2 
spectrum, the flow speed is ~0.1 km s~\ substantially faster 
than generally expected for ambipolar diffusion. If such a flow 
operates for 0.1-1 Myr, it may bring more mass into a dense 
core than does ambipolar diffusion, thus reducing the time of 
core condensation. Shorter core formation times would help 
account for the relatively low proportion of starless dense cores 
observed in Taurus and Orion, which suggest that starless cores 
are detectable for less than the few Myr expected from am- 
bipolar diffusion (Jijina, Myers, & Adams 1998). 

Turbulent cooling flows should also bear on the line 
width-size relations, in which nonthermal line widths increase 
with map size and decrease with mean map density (Larson 
1981; Myers & Goodman 1988; Falgarone, Puget, & Perault 
1992; Goodman et al. 1998). Since turbulent cooling flows 
increase the density of gas with low-velocity dispersion, they 
should tend to decrease the line width more on small size scales 
than on large size scales. If such flows are sufficiently prevalent, 
they should contribute to the trends evident in the line 
width-size relations. 
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On the other hand, much more work is needed to fully de- 
scribe the flows outlined here and to compare their properties 
with observations. Our treatment ignores the role of the static 
magnetic field in retarding inward motion and thus requires 
that the flow be primarily along field lines. It neglects the 
excitation of turbulence in the inner region by the flow and 
exaggerates the flow pressure gradient by considering a uniform 
inner region and equipartition turbulence in the outer region. 
We consider only one type of "turbulence," a linear superpos- 
ition of MHD waves, and only one mechanism of dissipation, 
ion-neutral friction, while wave steepening and nonlinear ef- 
fects should also be important. We ignore the role of self- 



gravity, while in reality turbulent cooling might precede and 
enhance gravitational infall. Finally, it is desirable to provide 
a space-time description of flow density and velocity to allow 
detailed comparison with observations. These topics are suit- 
able for future papers, and some are discussed by Lazarian & 
Myers (1998). 

We thank Fred Adams, Bruce Draine, Bruce Elmegreen, 
Charles Gammie, Alyssa Goodman, and Mario Tafalla for help- 
ful comments and Frank Shu for detailed comments and 
suggestions. 
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